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The invariant mass spectrum of the r]'-K + -K~ final state produced in two-photon collisions is ob- 
tained using a 673 fb" 1 data sample collected in the vicinity of the T(4S) resonance with the Belle 
detector at the KEKB asymmetric-energy e + e~ collider. We observe a clear signal of the rj c (lS) 
and measure its mass and width to be M(rj c (lS)) = (2982.7 ± 1.8(stat) ± 2.2(syst) ± 0.3{model)) 
MeV/c 2 and T(rj c (lS)) = (37.8t (stat) ± 2.8(syst) ± lA(model)) MeV/c 2 . The third error is an 
uncertainty due to possible interference between the r; c (lS) and a non-resonant component. We 
also report the first evidence for t;(1760) decay to n'ir + ir~; we find two solutions for its parameters, 
depending on the inclusion or not of the X(1835), whose existence is of marginal significance in our 
data. From a fit to the mass spectrum using coherent X(1835) and 77(1760) resonant amplitudes, 
we set a 90% confidence level upper limit on the product T 11 B{n'-K + -R~) for the X(1835). 

PACS numbers: 13.25.Gv, 14.40.Gx, 13.66Bc, 12.38.Qk 



I. INTRODUCTION 

As the lowest charmonium state, the r/ c (lS) meson 
plays an important role in tests of QCD. However, 
even its main parameters, such as the mass, width 
and two-photon width, have not been well measured 
and the measurements that have been reported show 
a large scatter of values [l|. Discrepancies among 
measurements for the t] c (1S) product of the two- 
photon width and decay branching fraction into four- 
meson final states were confirmed earlier 0]. A re- 
cent measurement of the r] c (lS) that found a signif- 
icant interference between the 77 C (1S) and the non- 
resonant background [1] may have clarified the rea- 
son for discrepancies among 77 C (1S) parameter mea- 
surements Q. Significant model-dependent uncer- 
tainty in the measurement of the rj c (lS) product 
branching fractions due to interference between the 
77c (IS) and a non-resonant component has also been 
studied in B — > Krj c (\S) decays 

The A(1835) resonance was observed and con- 
firmed recently by the BES collaboration in J/ip — > 
7X(1835) decays where A(1835) j/tt+ti- [|, with 
mass M = (1836.5 ± 3.0±^) MeV/c 2 and width 
r = (190 ± 9tU) MeV/c 2 . A variety of specula- 
tions on the nature of the A71835) have been re- 
ported, including baryonium [7J with sizable gluon 
content glueball [9hlll1|. and a radial excita- 
tion of the 77' [H EH. The BES experiment has 
suggested that the X(1835) may be related to the 



pp threshold enhancement seen in J/yj — > jpp de- 
cays [13, EH- An additional structure, the 77(1760), 
was observed in the radiative J/ip decays to jpp and 
jujoj by MARKIII and DM2 [H and to jujlj 
and 7777r+7r- by BES p3|. The 77(1760) state has 
been proposed as a mixture of a gluonic meson with 
a conventional qq state fliij ). rather than a pure qq 
meson, and this hypothesis is supported by a BES 
analysis of J/ip — > jlou) decays [18J. Hence, an in- 
vestigation of the nature of both the A(1835) and 
77(1760) is of interest [20|. In radiative J/ip decays, 
hadrons are produced via two gluons; thus, the pro- 
duction of final states with a gluon-enrichcd com- 
ponent is expected to be enhanced. In light of the 
similar structure of the two-photon and two-gluon 
couplings, a comparison of the 77 width of a meson 
to its production rate in radiative J/ip decays can 
provide information on its quark and gluon compo- 
sition. The two-photon coupling to the gluonic com- 
ponent is expected to be very weak so measurements 
of two-photon widths can help clarify the nature of 
the A(1835) and 77(1760). 

In this paper, we report the first observation of 
i]'tt + tt~ production in two-photon collisions using 
a 673 fb" 1 data sample (605 ftr 1 on the T(4S) 
resonance and 68 fb -1 at 60 MeV below the res- 
onance) accumulated with the Belle detector 1211 
at the KEKB asymmetric-energy e + e~ collider [22| . 
We measure parameters of the 77 C (1S), provide first 
evidence for 77(1760) — >• ifir + n~ decay, and give lim- 
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its on the two-photon production of the X (1835). 



II. DETECTOR AND MONTE CARLO 
SIMULATION 

The Belle detector is a largc-solid-angle magnetic 
spectrometer that consists of a silicon vertex detec- 
tor, a 50-layer central drift chamber (CDC), an array 
of aerogel threshold Cherenkov counters (ACC), a 
barrel-like arrangement of time-of-flight scintillation 
counters (TOF), and an electromagnetic calorimeter 
comprised of Csl (Tl) crystals (ECL). These detec- 
tors are located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. An iron 
flux return located outside the coil is instrumented 



to detect K\ mesons and to identify muons [21 1. 

Monte Carlo (MC) events of the two-photon pro- 
cess 7*7* — > ifir + ir~ are generated with the TREPS 
code [23| based on an Equivalent Photon Approxi- 
mation (EPA) [23, where the 77' decays generically 
according to the JETSET7.3 decay table An 
isotropic phase space distribution is assumed for 
r) c (lS), 77(1760) and X(1835) decays to the three- 
body 77'7r + 7r~ final state. The GEANT-based simu- 
lation package [26| with trigger conditions included 
is employed for the propagation of the generated par- 
ticles through the Belle detector. 



III. EVENT SELECTION 

The T) C (1S), 77(1760) and X(1835) (collectively de- 
noted as R) candidates are reconstructed from the 
decay chain R — > ifir + Tr~ , 77' — > r/TT + Tr~, and 77 — > 
77. Two photons and two tt + tt~ pairs are detected 
in the final state. 



A. Selection criteria 

At least two neutral clusters and four charged 
tracks with zero net charge are required in each 
event. Candidate photons are neutral clusters that 
have an energy deposit greater than 100 MeV in the 
ECL and are not near any of the charged tracks. 
The polar angle of the charged tracks, i.e., the angle 
with respect to the direction opposite the positron 
beam axis in the laboratory system, must satisfy 
cos6> G [-0.8660, +0.9563]. To enhance the detection 
efficiency for low momentum charged tracks, loose re- 
quirements on the impact parameters perpendicular 
to (dr) and along (dz) the beam line from the interac- 
tion point are applied: dr < 5 (< 3, < 2, < 1) cm and 



|g?;z|<5(<5, < 4, <3) cm for the track transverse 
momentum p t < 0.2 (e [0.2,0.3], € [0.3,0.4], > 0.4) 
GeV/c. The scalar sum of the absolute momenta for 
all the charged tracks and neutral clusters and the 
sum of the ECL cluster energies in the laboratory 
system are required to be p sum < 5.0 (< 5.5) GeV/c 
for the rj'ir + iT~ system in the mass region below 2.7 
GeV/c 2 (in the r) c (lS) region) and E sum < 4.5 GeV. 

Events with an identified kaon (K^ or K s — > 
7r + 7r~) or proton are vetoed. For charged tracks, 
information from the ACC, TOF and CDC is com- 
bined to form a likelihood C for hadron identifica- 
tion. A charged track with the likelihood ratio of 
£k/(£tt + £-k) > 0.8 is identified as a kaon; one 
with £„/(£„ + Ck) > 0.2 as a pion. With these 
loose requirements, the efficiency for pion identifi- 
cation is about 99%. A proton is identified by the 
requirement C P /(C P + Ck) > 0.95. K s candidates 
are reconstructed from a pair of charged pion tracks 
with invariant mass within 16 MeV/c 2 (3a) of the 
nominal K s mass. 



The 77 from 77' — > i]tt + tt~ decay is reconstructed via 
its two-photon decay mode, where the two-photon 
invariant mass is in the window M 11 G [0.524, 0.572] 
GeV/c 2 (±2tr of the nominal 77 mass). To suppress 
background photons from ir° decay, we exclude any 
photon that, in combination with another photon 
in the event, has an invariant mass within the win- 
dow |M 77 - m w o\ < 18 MeV/c 2 . The two-photon- 
energy asymmetry, A sym = \E^i~E 72 \/(E 7l + E 72 ), 
is required to be less than 0.8 to suppress the fake 
77 combinatorial background. The 77' candidate is 
reconstructed from the 77 candidate and the ir + n~ 
track pair that results in an invariant mass within 
M v „+„- G [0.951,0.963] GeV/c 2 (±2cr of the nomi- 
nal 77' mass). To improve the momentum resolution 
of the 77 and 77', a mass-constrained fit to the 77 and 
two separate fits to the 77' (one with a constrained 
vertex and the other with the mass constrained to 
the 77') arc applied. 



The rj'n + -K~ candidates are reconstructed by com- 
bining the 77' candidate and the remaining 7r + 7r~ 
track pair. For multi-candidate events, the candidate 
with the smallest \m n '°m the 77' mass-constrained fit 
is selected. For ij'tt^t:^ combinations with invariant 
mass W = 1.84 (2.98) GeV/c 2 , 19% (7.1%) of the 
signal MC events have more than one candidate per 
event, from which the correct candidate is selected 
98% (91%) of the time. 
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Background and optimization for | Pt I 
requirement 



Signal and non-resonant events can be produced 
in two-photon collisions via the processes e + e~ — > 
e + e~R and e + e~ — > e + e~ rf V + ir~ , respectively, 
where quasi-real photons are emitted from the beam 
e + and e~ particles at small angles with respect to 
the beam line. These events tend to carry small 
transverse momentum | ^ p t * | , which is determined 
by taking the absolute value of the vector sum of the 
transverse momenta of ?/ and the 7r + 7r~ tracks in the 
e + e~ center-of-mass system. 



The r/-sidcband, denoted rj'-sdb, arises from 
?77r + 7r~7r + 7r~ and 777r + 7r~7r + 7r~ (without 77) combi- 
nations that survive the rj selection criteria except 
that the rjir + Tt~ combination whose mass is nearest 
that of the vj lies between 0.914 and 0.934 GeV/c 2 
or between 0.98 and 1.0 GeV/c 2 . Similar events 
with an r/Trir mass within the rj' acceptance window 
form a featureless background denoted b\ in the R- 
candidatc sample. The rj'ir + Tr^ X background, de- 
noted &2 , has additional particles in the event beyond 
the R candidate. Other non-exclusive backgrounds, 
including those arising from initial state radiation, 
are found to be negligible. 

Significant background reduction is achieved by 
applying a \J2Pt\ requirement. The \J2Pt\ distri- 
bution for the signal peaks at small values, while that 
for both backgrounds decreases toward \J2,Pt\ = 
due to vanishing phase space [27 1. 



The r) c (lS) state is well established 0, H [1| and 
its signal yield in our data sample is large. We uti- 
lize a control sample of ?/7r + 7r~ candidates from half 
the data, with W between 2.6 and 3.4 GeV/c 2 , to 
establish the | Pt I requirement under the assump- 
tion that the | ^2 Pt distribution is similar for events 
with W < 2.2 GeV/c 2 . The rj'-sdb events from the 
full data sample are added to this control sample un- 
der the assumption that their ^2 Pt distribution is 
similar to that of the 61 background so that the sig- 
nal fraction in this control sample is close to that in 
the W mass region below 2.2 GeV/c 2 in the full data 
sample. We use the relative statistical error for the 
rj c (lS) yield in fitting the f/7r + 7r~ mass spectra to 
optimize the requirement. The requirement 

I J2Pt\ < 0-09 GeV/c (p t -balanced) is applied to the 
i?-candidate sample since it minimizes this relative 



IV. BACKGROUND ESTIMATION 

The 61 component in the r]'ir + TT~ mass and J^Ptl 
distributions are determined in the fits to the rj - 
sdb events (normalized) in the p t -balanced and pt- 
unbalanced (see below) samples, respectively. The 
residual bi component in the final i?-candidate sam- 
ple can be separated using the | ^2 Pt distribution. 
By doing so, its distribution in ?/7r + 7r~ mass is de- 
termined. FigureQ]shows the | J2Pt\ distribution for 
signal MC events and data in the mass region below 
2.2 GeV/c 2 . 




error. 



0.1 0.15 0.2 

lEpi l in GeV/c 

FIG. 1: The |X^Pt*| distributions for the mass region 
below 2.2 GeV/c 2 the data sample. The data points with 
error bars are from the r/'-K + -K~ -candidate sample before 
the I ^2Pt requirement, the thick-solid histogram is the 
best fit, the thin solid histogram is the signal component, 
the thin-dashed curve is the 61 component (whose shape 
is taken from the r/'-sdb sample), and the thin-dotted 
curve is the &2 component. 

A p t -unbalanced data subsample, in which the 
backgrounds dominate over the signal, is selected 
with the requirement \ J2Pt\ £ [0.15,0.2] GeV/c. 
The 7/7T7T mass distribution of this p t -unbalanced 
subsample is fit to two separate background func- 
tions, one for the 61 component with its yield and 
shape fixed at the values determined using the cor- 
responding r/'-sdb sample and the other for the &2 
component with its yield y un bai arid shape parame- 
ters allowed to float. We use the same shape for the 
62 component in the later fit to the rj'Tr + TT~ mass 
spectrum for the final i?-candidate sample. Here, 
the assumption of the same shape in the invariant 
mass distribution for the 62 component in the pt- 
balanccd and -unbalanced samples is implied. In the 
fit shown in Fig. [TJ the signal function for R and 
non-resonant events is defined by a histogram of the 
signal MC events with its shape parameters fixed 
but yield floated; the b\ component is described by 
a threshold function with its yield and shape pa- 
rameters fixed; the 62 component is described by 
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a quadratic function with its yield and shape pa- 
rameters floated. Here, the quadratic function for 
the b 2 is constrained to the origin, since b 2 back- 
ground events selected as r)'n + ir~~ with missing X 
should have non-zero transverse momentum. From 
the fit, we obtain the b 2 yields in the p 4 -balanced and 
-unbalanced subsamples; the ratio of these yields is 
Ubai/Uunbai = 0.723 ± 0.043. (The corresponding 
b 2 yield ratio for 2.6GcV/c 2 < W < 3.4GeV/c 2 
is 0.93 ± 0.11.) The b 2 yield y' bal in the t/tt+tt- 
mass spectrum for the final i?-candidate sample is 
obtained from the yield y' unba i scaled to this ratio. 

The invariant mass distributions for the r)'ir + 'K~ 
candidates, as well as those for the b\ and b 2 back- 
grounds, are shown in Fig. [5] In addition to the 
prominent r\ c (\S) signal, an enhanced shoulder is ev- 
ident in the mass region below 2 GeV/c 2 in the bl- 
and b 2 -subtracted histogram of Fig.[3Jb). The robust 
enhancement is also seen in the ?/7r + 7r~ yields ex- 
tracted from fitting the | ^2 P* I distributions in each 
sliced mass bin, shown as data points with error bars 
in Fig.gfb). 




1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 3.4 

M(h'tiV) [GeV/c 2 ] 




1.5 1.75 2 2.25 2.5 2.75 3 3.25 

M(h'tiV) [GeV/c 2 ] 

FIG. 2: Invariant mass distribution for the ti'tx + tx~ can- 
didates, (a) The open histogram represents the data; the 
horizontal (vertically) hatched histogram is the contri- 
bution from the 61 (62) background, (b) The histogram 
shows the data after subtraction of both the 61 and 62 
background components; the points with error bars are 
the r] / 7v + Tv~ yields extracted from fitting the | X]Pt*l dis- 
tribution in each sliced mass bin. 



V. FITTING MASS SPECTRUM 

The cross section of R production in the two- 
photon process e + e~ — > e + e~ R is approximated by 

<?(e + e- -> e+e-R) = J a 7 ^ R (W)^dW, (1) 

where the two-photon luminosity function ^27 is 
calculated in the EPA using TREPS and the cross 
section Ujj^ji(W) for C-even resonance production 
with zero spin is described by a Brcit-Wigncr (BW) 
function f B w{W) 

a^ R (W) = iW(w0-r 77 

(W 2 — M 2 ) 2 + Y 2 M 2 ' 1 ' 

where M, T and T 77 are the mass, total width and 
two-photon decay width of the R, respectively. 

The signal yield n s , M and T are extracted by 
maximizing the extended likelihood function, 

£ = W] Y[[n s -f s (u f ,M,r) 

i=l 

3 

+ nb ' k -/mK;pm)L ( 3 ) 

k=l 

where n s (nb,k) is the number of signal (fc-th back- 
ground component) events, N is the total number 
of candidate events, i is the event identifier and Ui 
is the measured invariant mass for the i-th event. 
The probability density function (PDF) f s for the R 
signal is a BW function convolved with mass reso- 
lution after corrections for -^pp- and the efficiency. 
The fc-th background's PDF and its parameters are 
denoted by /^j. and Pb,k, respectively. In the fit, 
n s , M and T for the signal are allowed to float un- 
less stated otherwise; and pb,k for non-resonant 
background {NR) are allowed to float while those 
for the bi and b 2 backgrounds are fixed. Two dis- 
tinct fits are performed: in the lower mass region 
1.4GcV/c 2 < W < 2.7GcV/c 2 where the NR (as 
well as bi and b 2 ) background component is de- 
scribed by a threshold function [2^] with a reasonable 
description of the threshold effect, and in the higher 
mass region 2.6GcV/c 2 < W < 3.4GeV/c 2 (near 
the r] c (IS 1 )) where all the background components 
arc described by an exponential of a third-order poly- 
nomial. 



6 



The evaluation of the significance of any marginal 
R signal in the lower-mass fit is sensitive to the as- 
sumed background shape. We have examined results 
of various fits with different descriptions of the back- 
ground: (1) one threshold function for a sum of all 
three background components (i.e., b\, b 2 and NR); 
(2) two separate threshold functions, one for b\ and 
the other for 62 plus NR; (3) three separate threshold 
functions, one each for 61, 62 and NR, respectively; 
(4-6) the three background functions defined above, 
in each case convolved with a mass resolution func- 
tion after corrections for the two-photon luminosity 
and efficiency We fit the i]'ir + Tr~ mass spectrum for 
a possible 7/(1760) signal in the mass region below 
2.7 GeV/c 2 using the six different background mod- 
els described above. Option (3) provides the small- 
est statistical significance for a signal resonance, and 
is conservatively chosen for the background descrip- 
tion. 

The product of the two-photon decay width and 
the r]'Tt + TT~ branching fraction for the R is deter- 
mined as: 



CM 

o 



> 




M(n 7t + 7t ) [GeV/c 2 ] 



FIG. 3: The invariant mass distribution for r]'ir + n~ 
candidates in the lower-mass region. The points with 
error bars are data. The thick solid line is the best fit; 
the thin solid line is the total background. The thick 
dashed line is the fitted signal for the 7/(1760) . The thin 
dashed, dot-dashed and dotted lines are the NR, bi and 
62 background components, respectively. 



r 77 /3(i? -> t/tt+tO = 



J f B w(W)^^le(W)dW 



(4) 



where the efficiency e includes the branching frac- 
tions for B(rf — > rjTr + ir~) and B(r/ — > 77). 



A. Results of the 7/(1760) fit 

We assume that only one resonance is produced 
in the mass range below 2.7 GeV/c 2 and that there 
is no interference between the signal and NR com- 
ponents. Figure [3] shows the results of the fit for 
the decay R — > ii'h + tt~. A signal with a yield 
n s = 465^24 and a statistical significance of 4.8a is 
found with mass M = (I768IJ5) MeV/c 2 and width 
r = (224rj 2 ) McV/c 2 ; we denote this as 7/(1760). 
The statistical significance, in units of standard de- 
viation (a), is calculated using the x 2 distribution 
—2-ln(£o/C m ax) with iVdof degrees of freedom. Here, 
£max and £0 denote the maximum likelihood with 
the signal yield floating and fixed at zero, respec- 
tively, and iVdof = 3 is the difference in the number of 
floating parameters between the nominal fit and the 
fit with the signal yield fixed at zero. The product 
of the two-photon decay width and branching frac- 
tion is determined to be r 77 /3(7/(1760) — > ?/7r + 7r~) = 
(28.21™) cV/c 2 . 



B. Results of the A(1835) fit 

According to existing observations [l|, |(| , two reso- 
nances, A(1835) and 7/(1760), have been reported in 
the lower mass region above the 7/ '7r + 7r~ threshold. 
Assuming that both X(1835) and r/(1760) have the 
same spin-parity of J PC = 0~ + , the effect of inter- 
ference between these two states must be taken into 
account in any attempt to extract a signal yield for 
the X(1835). Each resonance is described by a BW 
amplitude: 



gBw(W) 



(W 2 - M 2 ) + iTM ' 



(5) 



and the amplitude for the two interfering resonances 
is written as 

M{W) = A± ■ g BW i(W) + A 2 ■ g BW2 (W) ■ e # , (6) 

where <f> is the relative phase between the two reso- 
nances and A\ and A2 are normalization factors. 

Under the assumption of non-interference between 
the R and NR components, a fit with the A(1835) 
and 7/(1760) signals plus their interference is per- 
formed to the lower-mass events. Here, the X(1835) 
mass and width are fixed at the BES values @. 
We find two solutions with equally good fit qual- 
ity and the same 77(1760) mass and width; the re- 
sults are shown in Fig. 01 In cither solution, the 
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statistical significance is 2.9a for the X(1835) and 
4.1(7 for the 7/(1760). The relative phase between the 
two resonances is determined to be <j>\ = (287+5x)° 
for the constructive-interference solution and 4> 2 = 
(139lg 9 )° for the destructive-interference one. The 
signal yields for the two solutions are determined to 
be Yi = 332±^° and Y 2 = 632^^4 for the x(1835) 

and Yi = 52±|g and Y 2 = 3151/^ for the 7/(1760). 
The fitted mass and width of the 7/(1760) are M = 
(1703±if) McV/c 2 and T = (42^) McV/c 2 . Upper 
limits on the product r 77 £>(7/'7r + 7r~) for the X(1835) 
at the 90% confidence level are determined to be 
35.6 eV/c 2 and 83 eV/c 2 for the constructive- and 
destructive-interference solutions, respectively. The 
upper limit for the signal yield at 90% confidence 
level is determined by integrating the likelihood dis- 
tribution convolved with a Gaussian function to in- 
clude the systematic error. 

Another fit without interference between the res- 
onances is performed to examine the significance of 
the X(1835) signal. The statistical significance from 
the fit with two incoherent resonances is found to 
be 3.2cr for the X(1835) and 4.4cr for the r/(1760). 
The 7/(1760) mass and width are fitted to be M = 
(1707.7±7.o) MeV/c 2 and T = (45^) MeV/c 2 , re- 
spectively. The products of the two-photon decay 
width and the branching fraction for the X(1835) 
and 7/(1760) decays to 7/'7r + 7r~ are estimated as 
r 77 /3(X(1835) -> 7/7T+7T-) = (23.1±e.6) eV/c 2 and 
r 77 £(r/(1760) -> rt'n+n-) = (6.7t 2 J) cV/c 2 . The 
inclusion of the interference only mildly improves 
the fit. The statistical significance of the interfer- 
ence term, defined as \J—2 ln(£ no /£ ycs ), is 0.69cr, 
where £ yos (£ no ) is the likelihood value of the fit 
with (without) interference. There is a minor differ- 
ence in the 7/(1760) mass and width between the two 
fits with and without interference. The statistical 
significance of the 7/(1760) mass difference between 
the fit result and the world-average value [l[ is calcu- 
lated as y— 2 ln(£fi xe d/£free), and is found to be 2.6<r 
(3.1cr) for the two coherent (incoherent) resonances. 
Here, £fi X ed and £f ree are the likelihood values of the 
fits with the 7/(1760) mass fixed at the world-average 
value and floating, respectively. 

In the determination of the X(1835) and 7/(1760) 
significances, we have examined the effect of uncer- 
tainties of the following factors: (1) the X(1835) 
mass or width varied by ±ler; (2) a background 
fluctuation by changing the fit region; (3) a back- 
ground function that uses three threshold functions 
convolved with two-photon luminosity, efficiency and 
mass resolution; (4) a fluctuation in the b\ compo- 
nent by moving the rj'-sdb selection mass window; (5) 
a variation of ±ler in each of the background function 
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FIG. 4: Results of a combined fit for the X(1835) and 
7/(1760) with interference between them. The points with 
error bars are data. The thick solid line is the fit; the thin 
solid line is the total background. The thick dashed (dot- 
dashed, dotted) line is the fitted signal for the 7/(1760) 
(X(1835), the interference term between them). The 
thin dashed, dot-dashed and dotted lines are the NR, b± 
and 62 background components, respectively. The upper 
(lower) panel represents the solution with constructive 
(destructive) interference. 



parameters for the b\ or b 2 components. The fits of 
two incoherent resonances are performed under the 
variations listed above. The lowest (highest) signifi- 
cance 3.9cr (5.0a) for the 7/(1760) is obtained with the 
X(1835) width increased (decreased) by la, while 
the significances under the rest of variations are com- 
patible with the values from the incoherent fit of 3.2(7 
for the X(1835) and 4.4cr for the 7/(1760). To ensure 
reliable estimation for the X(1835), a fit with float- 
ing masses and widths for both the X(1835) and the 
7/(1760) is performed. The yields, masses and widths 
are fitted to be Y = (444 ± 158), M = (1833 ± 30) 
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MeV/c 2 and T = (202±66) MeV/c 2 for the Jf (1835) 
and Y = (104 ± 75), M = (1706.9 ± 8.3) MeV/c 2 
and r = (40 ± 36) McV/c 2 for the 77 (1760). In all 
variations, the fitted parameters for the A"(1835) are 
consistent with those in the BES experiment. 



C. Angular distribution 




0.2 0.4 0.6 0.8 t 1 

|COS(0*)| 

FIG. 5: Angular distributions in the 77 rest frame for 
two mass regions. The solid circles are for the X(1835) 
and 77(1760) region; the open circles are for the N R com- 
ponent in the upper sideband region. The yield in each 
bin is corrected for the efficiency and normalized to the 
sum of the corrected yield. The horizontal line at y = 0.1 
represents an isotropic MC distribution. 



of pseudoscalar quantum numbers for the »7(1760) 
and X(1835). However, a possible non-flat distribu- 
tion for the NR will influence the distribution for 
the R component; thus, a plausible J p value for 
each R should be examined with the NR, compo- 
nent subtracted once the existence of the 77(1 760) 
and AT(1835) production is clarified. No significant 
intermediate state is seen in the mass region below 
2.04 GeV/c 2 . However, a minor contribution from 
another J p = 2~ resonance [3(| cannot be ruled 
out. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Invariant mass of n*n ' 



FIG. 6: Invariant mass distribution of ir tt for the 
rj'ir + Tv~ events selected in the upper sideband region of 
2.2 GeV/c 2 < W < 2.7 GeV/c 2 , where a large signal for 
/2(1270) — > tv + tt~ decays is shown. 



We examined the distribution of 9* , the angle be- 
tween the 77' momentum and the beam direction in 
the 77 rest frame. The angular distribution is de- 
termined from R and NR yields extracted from fit- 
ting the I ^2 Pt I distribution sliced into ten angular 
bins for the mass region of the AT(1835) and 77(1760) 
(W < 2.04 GeV/c 2 ) and its upper sideband (W € 
(2.2,2.7) GeV/c 2 ). The distribution in the upper 
sideband region shows forward and backward peaks 
characteristic of a higher-angular-momentum com- 
ponent, which indicates strong contributions from 
the rj'f 2(1270) production (see Fig. Indeed, a 
large /2(1270) signal is observed in the 7r + 7r~ in- 
variant mass distribution for the r)'n + TT~ events se- 
lected in that region, as shown in Fig. [5] The domi- 
nant ?//2(1270) component in the upper sideband re- 
gion shows interesting dynamics with a broad struc- 
ture with favored quantum numbers J p = 2 + . A 
nearly isotropic distribution in the mass region be- 
low 2.04 GeV/c 2 after the efficiency correction (with 
X 2 /Ndof = 9.9/9) is compatible with the assumption 



D. Results of the r) (lS) fit 

We first assume that there is no interference be- 
tween the ?] C (1S) and the NR background. Fig- 
ure [7] shows the ti't: + t:~ invariant mass distribu- 
tion for the candidates with mass greater than 2.6 
GeV/c 2 together with the fitted signal and back- 
ground curves. The ?7c(15) mass and width are de- 
termined to be M = (2982.7 ± 1.8) MeV/c 2 and 
r = (37.8±|;|) McV/c 2 . The product of the two- 
photon decay width and branching fraction for the 
77c (15) is calculated using Eq. (HJ). Using the fitted 
77c (15) signal yield of n s = 486lgg, we determine 
r 77 /3(77 c (15) ?7'7r+7r-) = (50.5±ti) eV/c 2 . 

We now address the effect of possible interference 
between the ?7c(15) resonance, hereafter referred to 
as i?, and the non-resonant component. A precise de- 
scription of the data in this case is impossible with- 
out a good understanding of the background. As 
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FIG. 7: The invariant mass distribution for the rj'n + n 
candidates in the mass range above 2.6 GeV/c 2 . The 
points with error bars are data. The thick-solid line is 
the fit; the thin-solid line is the total background. The 
thick dot-dashed line is the fitted signal for the r] c (lS). 
The thin dashed, dot-dashed and dotted lines are the 
NR, bi and 62 background components, respectively. 

discussed in section V-C, the NR component in the 
mass region above 2.2 GeV/c 2 has a contamination of 
events from non-0 _ production via two-photon pro- 
cesses. Although contamination is evident even in 
the rj c (IS) mass region, our data sample is insuffi- 
cient to determine the type and rate of production of 
the non-0~ states in this mass region. The NR com- 
ponent in our analysis can be subdivided into two 
types: one for the non-resonant final state (denoted 
as NR1) that interferes with R and the other for 
production of various non-0~ states (denoted NR2) 
that do not interfere with the R. The amplitude for 
R production with interference with the NR1 term 
is written as 

M(W) = A ■ g B w{W) ■ + A NR1 ■ g Nm (W), (7) 

where gsw is the BW function in Eq. (JSJ) , gNRi 
is assumed to be a real function for NR1, <f> is the 
interference phase, and A and Anri are normaliza- 
tion factors. Assuming that NR1 and A^i?2 have the 
same shape, the fitting function in Eq. ([3]) for the R 
and NR components — where R interferes with NR1 
but not with NR2 — can be expressed as 

/ = n s ■ f 8 (u; M, T) + n NR ■ /nr{w, pnr) + fint, (8) 

where the interference term is 

fint = 2 a/ (XNR ■ UNR '• fNR,{u] Pnr) 

■ y/n a -f s (u;M,T) ■ cos(9 + fa (9) 



with aNR = ni/riNR, unr. = n\ + ?i2, and n\ and 
ri2 are the number of NR1 and NR2 events, respec- 
tively. An intrinsic phase 9 is determined by the 
R mass, width and W value. The function /, in- 
cluding the fi n t term, is convolved with a mass res- 
olution function after corrections for dL^/dW and 
efficiency. The f s and $nr PDFs are normalized; 
the function fi nt is fully determined by the fit pa- 
rameters. 

To investigate the possible effect of interference 
with the NR component, a fit to the r] c (lS) sig- 
nal with interference with A^i?l but without inter- 
ference with NR2 is performed for various initial in- 
put values for the aNR and 4> parameters. For the 
t] c (1S), the fit gives two solutions with almost the 
same maximum likelihood value; the mass and width 
of the 7 lc (lS) are M = 2982.7 (2983.0) McV/c 2 and 
r = 36.4 MeV/c 2 at a NR = 0.01% (100%); these 
are quite consistent with the result of the fit with- 
out interference. The differences in the rj c (lS) mass 
and width with and without interference, AM = 0.3 
McV/c 2 and AT = 1.4 McV/c 2 , are taken as model- 
dependent uncertainties in the determination of the 
mass and width. However, the fits give very different 
values for the t] c (1S) yield. If, for example, aNR. is 
fixed at 100% in the fit with interference, the yields 
obtained are Y 1 = 854 ± 59 with fa = (-92 ± 5)° 
for destructive interference and Y2 = 264 ± 22 with 
fa = (91 ± 8)° for constructive interference, while 
the T)JlS) yield of the incoherent fit is 486±3g. A 
strong correlation between olnr and <j) 1S observed 
from the fits: fa and fa are close to 180° and —180° 
(90° and —90°), respectively, if aNR is close to zero 
(100%) . The insensitivity of the maximum likelihood 
value for the fits in the full aNR. region between zero 
and 100% and a strong correlation between aNR and 
<j) imply large uncertainties in the determination of 
aNR, <P an d t ne strength of the interference term. 
With an additional error source from the interfer- 
ence term, the r) c (lS) yield has also a large uncer- 
tainty ranging from 264 ±22 to 486 ±40 for construc- 
tive interference and from 486 ± 40 to 854 ± 59 for 
destructive interference depending on the true aNR 
and 4> values. Our fit results, as well as the absence 
of any visual asymmetry in the rj c (lS) line shape 
in the data, indicate that the interference term can- 
not be determined without independent information 
on the NR1 component such as its angular distribu- 
tion in the ri c (lS) sideband mass region. The mea- 
sured mass and width of the ri c (lS) have a marginal 
dependence on the interference, while the yield is 
strongly correlated with the interference component 
and, thus, cannot be determined precisely with the 
existing data sample. The situation would improve 
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if the interference effect were determined indepen- 
dently with a much larger data sample. 



VI. SYSTEMATIC ERRORS 

To examine a possible bias in the mass measure- 
ment for the decay R — > rj'ir + Tr~, a data sample of 
D° — > rj'Kg decays with K s —> ir + n~ is selected 
with tight mass window requirements for the rj and 
rf . The D° mass resulting from a fit of the invari- 
ant mass spectrum of ij'Kg is lower than its nominal 
value by 1.4 MeV/c 2 , which is taken as an uncer- 
tainty of the mass scale after a linear correction for 
mass value. The uncertainty in the width determina- 
tion can arise from a difference in the mass resolution 
between data and MC simulation. This is estimated 
by changing the mass resolution by ±1 MeV/c 2 and 
is found to be 2.0 MeV/c 2 for the r] c (lS) and 10 
MeV/c 2 for the 77(1760). Systematic errors on the 
mass, width and T J7 B product due to uncertainties 
in the NR background estimation are determined by 
varying the fit mass interval and | ^ jj t * requirement 
separately. The error contributions from uncertain- 
ties in determination of the b\ and hi backgrounds 
are minor for the ri c (lS) but are sizable in the mass 
region below 2 GeV/c 2 . The uncertainties in the res- 
onance parameters, estimated by varying the shape 
parameters and yields of the b\ and 62 backgrounds 
by ±ler and added in quadrature, are taken as the 
corresponding errors for the A(1835) and 77(1760) , 
respectively. 

There are additional sources of systematic errors 
in the r 77 £> product determination. The trigger ef- 
ficiency for four-track events is relatively high be- 
cause of redundant two-track and multi-track trig- 
gers in the Belle first-level trigger. From the trigger 
simulation program, the difference in the efficiency 
with and without both trigger conditions satisfied 
is found to be 1% (2.7%) at an invariant mass of 
2.98 (1.84) GeV/c 2 ; this is included as a systematic 
error. The efficiency for the pion identification, de- 
termined by using the inclusive D* sample, is lower 
than that from MC simulation by (1.40 ± 0.64)% 
for the T) C (1S) and (0.02 ± 0.60)% for the 77(1760), 
and the corresponding contributions to the system- 
atic error are 1.5% and 0.6%, respectively. The re- 
construction efficiency for 77 — > 77 is studied with an 
inclusive 77 sample, and its deviation from the MC 
simulation plus its error in quadrature is 4.9%. The 
uncertainty in the track reconstruction efficiency is 
5.5% and that of the 7r°-veto requirement is 3%. The 
accuracy of the two-photon luminosity function cal- 
culated by the TREPS generator is estimated to be 



about 5% including the error from neglecting radia- 
tive corrections (2%), the uncertainty from the form 
factor effect (2%), and the error of the total inte- 
grated luminosity (1.4%) [23|. The background con- 
tribution from the initial-state radiation processes is 
negligible 0] . Furthermore, the run-dependent back- 
ground conditions add an additional uncertainty of 
3% in the yield determination. A dominant source 
of systematic errors for the A(1835) yield is the un- 
certainty of its decay width. It is estimated to be 
18% by changing the width by ±lerr in the fit for 
the yield extraction. 

The systematic errors in the measurements of the 
mass and width for the r) c (lS) and 77(1760), as well 
as of the product F 77 B for each resonance, are sum- 
marized in Table |TJ 



VII. RESULTS AND DISCUSSION 

The results for the yields, masses and widths, as 
well as the product decay widths are summarized 
in Table HI for the r? c (15) and in Table HH] for the 
77(1760) and A(1835). 



The rj c (lS) mass and width are measured to be 
M = (2982.7 ± l.8(stat) ± 2.2{syst) ± 0.3(model)) 
MeV/c 2 and T = (37 .8±H(stat) ± 2.8(syst) ± 
1 A(model)) MeV/c 2 , and are consistent with the re- 
cent results from BES and Belle [Bj]. If we as- 
sume that there is no interference, the directly mea- 
sured product for the t] c (1S) decay width to t(tx + t;~ 
is determined to be r i7 B(rj c (lS) — > 7]'tt + tt~) — 
(50.5^4 2 ± 5.6) eV/c 2 , which is marginally consis- 
tent with the existing value (194 ± 97) eV/c 2 from 
the indirect measurements [l|. Instead of a direct 
reference to the world-average value for r 77 (r7 c (15)), 
we determine it from the ratio of T-y 7 T(K Ktt) /Ftotai 
= (0.407 ±0.027) kcV/c 2 to T(K 'Kit) /Ttotai = (7.0 ± 
1.2) x 10" 2 0, and obtain the width T 77 (t7 c (15)) 
= (5.8 ± 1.1) keV/c 2 with a smaller relative error. 
With that as an input, the branching fraction is cal- 
culated to be B(rj c {lS) rfir+%-) = (0.87±0.20)%, 
where both statistical and systematic errors are in- 
cluded. 

We report the first evidence for 77(1760) decay to 
77'7r + 7r~ and find two solutions for its parameters, 
depending on the inclusion or not of the A (1835), 
whose existence is marginal in our fits. The decay 
77(1760) — > rj'ir + ir~ is found with a significance of 
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TABLE I: Summary of systematic uncertainty contributions to the mass and width for the rj c (lS) and 77(1760) and 
to T 77 B for the r/ c (lS), 77(1760) and X(1835). 1-7? and 2-R denote one and two resonances in the fit, respectively. 



Source 




77(1760) 


A' (1835) 




1-R fit 






2-R fit 




A(M) (MeV/c 








Mass scale 


2.2 




1.3 




Background shape 


0.1 


8 


0.5 




rf sideband and b any 


0.0 


3.9 


0.2 




~)2Pt\ requirement 


0.4 


4.5 


0.6 




X(1835) Width 






0.9 




Total 


2.2 


10 


1.8 






A(r) (MeV/c 


*) 






Mass resolution 


2.0 


10 


1.5 




Background shape 


1.9 


7 


6 




77' sideband and b any 


0.02 


17 


7.1 




^2Pt\ requirement 


0.4 


14 


9 




X(1835) Width 






8 




Total 


2.8 


25 


15 




A(r 77 /3)/(r 77 B) (%) 


X(1835) Width 






16 


18 


Background shape 


4.6 


2 


13 


2.6 


rf sideband and b any 


0.03 


7.3 


15 


3.8 


Pt I requirement 


2.2 


0.6 


6.9 


6.3 


Trigger efficiency 


1 




2.7 




n ID efficiency 


1.5 




0.6 




77 rec. efficiency 






4.9 




Track rec. efficiency 






5.5 




7T° veto 






3 




Two-photon Luminosity 




5 




Run dependence 






3 




Total 


11 


13 


28 


22 



TABLE II: Summary of the results for the r] c (lS): M and 
r are the mass and width; Y is the yield; B is the branch- 
ing fraction for rj c (lS) — > rj'-K + Ti^ ; r 77 /3 is the product of 
the two-photon decay width and the branching fraction. 
The world-average values are shown for comparison. 



Parameters This work PDG 

Y 486t^ ± 53 

M, MeV/c 2 2982.7 ± 1.8 ± 2.2 2980.3 ± 1.2 

T, MeV/c 2 37.81:5 3 ± 2 -8 26.7 ± 3 

r 77 /3, eV/c 2 50.5lt 2 ± 5.6 194 ± 97 

£>, % 0.87 ± 0.20 2.7 ±1.1 



4.7ct including systematic error, with the assumption 
that the X(1835) is not produced; the 7/(1760) mass 
and width arc determined to be M = (1768±25 ± 10 ) 
MeV/c 2 and Y = (224±^ ± 25) MeV/c 2 . The 
fitted ?7(1760) mass is consistent with the existing 
measurements [l?], [H|. The product of the two- 
photon decay width and the branching fraction for 



the 77(1760) decay to ifn+ir is determined to be 
r 77 /3(77(1760) -> rfjT+TT-) = (28.2±£s ± 3.7) cV/c 2 . 
When the mass spectrum is fitted with two coherent 
resonances, the 77(1760) and X(1835), the 7/(1760) 
mass and width are found to be M = (1703 1 \l ± 1.8) 
MeV/c 2 and T = (42 j;^ ±15) McV /c 2 , and the sig- 
nal significances including the systematic error es- 
timated to be 4.1cr for the 77(1760) and 2.8<r for 
the X(1835). Upper limits on the product r 77 /3 
for the X(1835) decay to rj'Tr + ir~ at the 90% con- 
fidence level for two fit solutions are determined: 
r 77 /3(X(1835) -4 ?7'7r+7T-) < 35.6 eV/c 2 with 
4>i = (287^5 2 ) for constructive interference and 
r 77 /3(X(1835) t/'tt+tt-) < 83 eV/c 2 with 
4>2 = (139lg 9 ) for destructive interference. 

In summary, we report the first observation of 
rj'ir + Ti~ production in two-photon collisions. Wc 
measure the mass, width and the product of the two- 
photon width and the branching fraction for the 77 c . 
We also report the first evidence for the rj'iT + ii~ de- 
cay mode of the 77(1760). No strong evidence for the 
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TABLE III: Summary of the results for r?(1760) and X(1835): M and F are the mass and width; Y is the yield; 
r 77 B is the product of the two-photon decay width and branching fraction; Ygo and (r 77 0) 9O are the upper limits 
at 90% confidence level with systematic error included. The t;(1760) mass and width from the two-resonance fit 
with interference, as well as world average values, are shown for comparison. S is the signal significance including 
systematic errors. 



Parameter One resonance Two interfering resonances 



Reference 



Solution I 



Solution II 



M, MeV/c 2 
T, MeV/c 2 
Y 
Y 90 

r 77 /3, eV/c 2 
(r 77 S) 90 eV/c 2 
S, a 



X(1835) 



1836.5 (fixed) 
190 (fixed) 
± 73 632 



1836.5 ±3.0 



oon+140 
oo^_ 122 

< 650 
18.2l££±4.0 

< 35.6 

2.8 



Si ± 139 
< 1490 



190 ± 9±H [6] 



35 



+ 12 
-13 



±8 



< 83 



M, MeV/c 2 
T, MeV/c 2 
Y 

r 77 S, eV/c 2 
S, a 



1768^25 
224±« 2 ± 25 
4651^4 ± 60 

1-7.9 



77(1760) 

i™o4 



1703tii ± !-8 
42± 2 l 2 i ± 15 
52±^ ± 15 315± 22 | ± 8* 
3.0l 2 -2±0.8 18l^±5 
4.7 4.1 
(287lg)° (139l^ 9 )° 



1756 ± 9 [1] 
96 ±70 [1] 



X(1835) is found. 
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